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Abstract. The discrete-time quantum walk is a quantum counterpart of the random walk. It is expected that the 
model plays important roles in the quantum field. In the quantum information theory, entanglement is a key resource. 
We use the von Neumann entropy to measure the entanglement between the coin and the particle's position of the 
quantum walks. Also we deal with the Shannon entropy which is an important quantity in the information theory. 
In this paper, we show limits of the von Neumann entropy and the Shannon entropy of the quantum walks on the 
one dimensional lattice starting from the origin defined by arbitrary coin and initial state. In order to derive these 
limits, we use the path counting method which is a combinatorial method for computing probability amplitude. |||||||| 

1 Introduction 

The discrete-time (or coined) quantum walk (QW) has been extensively studied by many authors as a 
quantum analogue of the random walk [2l[6l[T5l[T6l[T9l[26]. As the random walk plays important roles in 
various fields, it is expected that the QW also plays such roles in the quantum field. In fact, benefits of 
using QWs have been shown in various applications, for example, quantum speed-up algorithm pilSl I 14112^ 
and universal quantum computation [T51[^ . 

In this paper, we give asymptotic behaviors of the von Neumann entropy and the Shannon entropy of 
discrete-time QWs on Z, where Z is the set of integers. Because entanglement does not appear in classical 
systems, it is an important concept for quantum information processing (see Nielsen and Chuang [23)). We 
use the von Neumann entropy of the reduced density matrix of the coin to quantify the entanglement between 
the coin and the particle's position. The Shannon entropy is one of a basic quantity in the information theory 
and it clarifies information included in the system. In the present paper, we calculate the entropies of the 
QWs starting from the origin with arbitrary coin and initial state from a path counting approach which 
is a combinatorial method for computing probability amplitude. By numerical simulations, Carneiro et 
al. [TT] studied the long-time asymptotic coin-position entanglement (aCPE) of QWs on various graphs, 
for examples, 1,1?, triangular lattices, cycles. Venegas-Andraca and Bose [5T also investigated the von 
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Neumann entropy on Z numerically. Using Fourier analysis techniques, aCPE of the Hadamard walk on Z 
for both localized (i.e., our setting) and non-localized initial conditions was analytically computed by Abal 
et al. [T]. In the similar technique, Annabestani et al. [8] gave an exact characterization of the aCPE of 
QWs on 1? . The evolution is determined by the tensor product of two one-qubit Hadamard operations. 
Liu and Petulante [3T] presented limit theorems for the von Neumann entropy of QWs on the A^-cycle. 
Bracken et al. [TU] numerically computed the Shannon entropy of the QW on Z defined by a coin which is a 
generalization of the Hadamard coin. The numerical result of Chandrashekar et al. |12j suggests that several 
properties of the Shannon entropy the QW on Z given by another generalization of the Hadamard coin. 

The rest of this paper is organized as follows. In Sect. 2, we give the definition of the QW. Results on 
the von Neumann entropy are presented in Sect. 3. Section 4 is devoted to the proof of our main result 
f Theorem 13. ip . In the final section, a long-time asymptotic for the Shannon entropy of the QW is shown. 



2 Definition 

The discrete-time QW is a quantum counterpart of the classical random walk with additional degree of 
freedom called chirality. The chirality takes values left and right, and it means the direction of the motion 
of the walker. At each time step, if the walker has the left chirality, it moves one step to the left, and if it 
has the right chirality, it moves one step to the right. Let define 



\L) = 



\R) = 



where L and i? refer to the left and right chirality state, respectively. 

Let U(2) denote the set of 2 x 2 unitary matrices. The time evolution of the QW on Z is determined by 



V 



a b 
c d 



e U(2), 



with a, 6, c, d G C and C is the set of complex numbers. The unitarity of U gives 
lap + = |c|2 + = 1, ac + 6rf = 0, c = -Ab, d = Aa, 
where z is the complex conjugate of z € C and A = det U = ad — be with | A| = 1. In particular, we write 



cost/ 
sin 9 



sin 9 

- cos 9 



where < 9 < Tr/2. When 9 = n/A, the QW is called the Hadamard walk. 

In order to define the dynamics of the model, we divide U into two matrices: 



P = 



with U ^ P + Q. The matrix P (resp. Q) represents that the walker moves to the left (resp. right) at each 
time step. Let S„(^, m) denote the sum of all paths starting from the origin in the trajectory consisting of I 
steps left and m steps right. In fact, for time n = I + m and position x = —I + m, we have 



a b 




■ " 





, Q = 


c d 



S„(/,to) 



P^^Q^'^P^^Q""^ ■ ■ ■ pi^-iQ^^-^-ipinQyyi„ 



where the summation is taken over all integers lj,mj > satisfying li + 
m, Ij + ruj = 1. We should note that the definition gives 

En+i{l,m) = P E„{1 - 1, m) + Q S„(/, m - 1) 

For example, in the case of Z = 3, m = 1, we have 

S4(3, 1) = QP^ + PQP^ + P^QP + P'^Q. 



+ h 



I, mi + ■ ■ ■ + m„ 
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The set of initial qubit states at the origin for the QW is given by 

$ = {(p = a\L) + I5\R) e C : lap + = 1} , 

where T is the transposed operator. The probabihty that a quantum walker is in position x at time n 
starting from the origin with </? G $ is defined by 

P(X„ = a;) = ||S„(/,TO)^||2, 

where n = I + m and x = —I + m. The probability amplitude ^n{x) at position x and time n is given by 

^n{x) = ^^{x)\L) + ^^{x)\R) = S„(Z, m)ip. (2.1) 

So p{x^ = x) = \^i:{x)\^ + M{x)\\ 

Let \^n) = • .,*„(-!), *„(0), *„(!), . . .]. The density operator at time n is given by /?„ = |*„)(^'„|. 
Entanglement for pure states can be quantified by the von Neumann entropy of the reduced density operator 

= Tip (pn), where the partial trace is taken over position. The associated von Neumann entropy at time 
n is 

S^ = -Tr {pllog.ip^J}. 

This is known as entropy of entanglement, which quantifies the quantum correlations present in the pure 
state. Let 

n n 

x=—n x=—n 

n 

x=—n 

Note that + ||^^||^ = 1 for any time n. The entropy of entanglement can be obtained after diago- 

nalization of p^. This is represented by the following 2x2 Hermitian matrix 

Let 

A„(p) = det{pl) = - K*^,*^)r 

Then the reduced entropy is expressed as 

Sn=- l0g2(?'n, + ) + r„ _ log2(r„ _)) , 

where r„,± are eigenvalues of p!^ and given by 



i±^/T^^iAj;p) 

rn,± = ^ . 

Here 01og2 is interpreted as 0. The determinant An{p) € [0, 1/4] quantifies the coin-position entanglement 
in the QW. The greater the value of A„(p), the greater the entanglement. If A„(/9) = 0, i.e., a product state, 
then = 0. If An{p) = 1/4, i.e., a maximally coin state, then = 1. We will compute the asymptotic 
value of as time n — >■ oo for the QW determined by J7 e U(2) with initial qubit state (fi G ^. 

On the other hand, the entropy of the reduced density matrix of the position at time n, = Ttc (Pn), 
also quantifies the entanglement between the coin and the walker's position. This is represented by oo x oo 
Hermitian matrix whose {x,y) element, p^(a;,y), is 



pP^ix,y) = ^ii{x)^^{y) + ^^{x)^S{y), 
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where a;, y € Z. The diagonal element becomes 

fP^{x,x) = + = P(X„ = x). 

Note that if \x\ > n or \y\ > n, then f)^{x^ y) = 0. Therefore the reduced entropy, 5^, can be expressed as 

n 

Sn = - ^ rn{x) logs {rn{x)) , 
x=—n 

where r„(a;) are eigenvalues of (2n + 1) x (2n + 1) Hermitian matrix with the element pP^{x,y) for x, y £ 
{— n, —n + 1, . . . , n}. We should remark that = for any n > 0, so we focus on S^. 



3 Results on entropy of entanglement 

In this section we present the following main result of aCPE for the QW on Z. Put 

11*^11= lim 11*^11, ||vl/«||= lim 11*^11, {^^,^^)= lim(vl/« vE-^). 



Theorem 3.1 When the QW is determined by U with abed ^ 0, we have 

ll^cx,ll 2 ' ' ^ V 2 2(1 + |6|)' 

w/iere F = aabl3 + aabji. Then 

lim S^ = - {roo,+ log2(roo, + ) + Too - log2(roo,-)} , 
n— ^oo 

where 



1± Vl-4Aoo(p) 
roo,± = ^ . 

and 

Aoo(p) = ii*^ini*sir-K*^,*S)r 

In particular, if ?7 = U{6), then the following result can be obtained. 

Corollary 3.2 When the QW is determined by U{9) with Q <0 < 7r/2, we have 

ii-rf 1,9 f. sin6'\ , ,9 sin6',^,, sin^fl — sin^) - — . „^ 
11*^11'= + — + 2cos^ (3.2) 

.^..2 sin0 2 /-I sin0\ 2 sin6'(l — siii^?) ^ — 



liuR ,T.Ls _ sing(l-sing) [|„,|2 1^,2 , sing (q^ + 5/3) + - q/3 | 
^ °°' 2COS0 V ' ^ cosg r 
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Reversible cellular automata were considered in 20 by the Fourier analysis. As a special case, the automaton 
includes the QW defined by U{9). Eqs. p.2p and p.3l) are equivalent to equations given in pp. 416-417 of [20]. 
The correspondence is H^"^!! o ||*^|| ||^'^(n)||, a O ai, j3 O Ur- Here we use a different 

approach based on a path counting. Abal et al. [1 investigated the case of U{'k/A) (i.e., Hadamard walk) 
for the localized initial condition (i.e., our setting) and the non-localized initial condition in the position 
space spanned by | ± 1). Noting that their notation ri.2 corresponds to our notation roo.±, their result for 
the localized case of = -^[cos^i, e'^^ sin6'i] g $ (stated at Appendix in [1]) is consistent with Corollary 
13.21 As for the QW on the line starting at the origin, Carneiro et al. |TT] studied the QW defined by 
U with a = ~d — ^fp, b — c = \J\ — p {p — 0.5 case is the Hadamard walk) and the initial condition 
If — ^[cos6'i, e*^^ sin 6*1] G They numerically found that the entanglement oscillates around an asymptotic 
value and the rate of convergence depends on the symmetry in the distribution of the QW. More symmetric 
distributions about the origin converge faster than asymmetric ones. Their limiting value for the Hadamard 
case is consistent with Corollarv 13.21 Recently, Alles et al. [3 studied an entanglement among two walkers 
on the line. We note that the model can be exactly mapped to our model by the following correspondence 
in their notations: 






^(Icft move) 



right move) 



4 Proof of Theorem 13.11 



In this section we assume abed ^ 0. We consider the following four matrices: 



P 



a b 




, Q 





c d 



R 



c d 




s 





a b 



Put X Ay ^ minja::, y}. For I A m > 1, we have 



I Am 



7=1 



2 \ 7 



/— l^^m — 1"^ ^Z — 7 ^ m — 7^ 



7— l/\7— 1 



07 



Aa7 



A6 



R+-S], 



by the path counting method [T71[TS]. Therefore, by noting Eq. (|2.ip . we have 

l^ra . ,,,2^7 n-l\fm-l 



^'^(x) = a'a^A"^ 

7=1 

/Am 

7=1 



l-l 



A =B 

a-f Ab 



771 — 7 ^ 1 ^ 

Aa7 b 



where A = aa + b(3, B — ca + df3. Let [x] denote the integer part of a; G M. We first consider Hvl/^ll case. 
For 1 < ? < [n/2], 



7=1 5=1 ^ 

P\A\^ I (\A\ 



\b\' 



7+5 



l-l\/'l-l\/'n-l-l 
j-l)[s-l)[ 7-1 



n-l-1 
S-1 



7^|ap 7 I |a| 



r rr^ + Q -7 rr^ + e + r:^ + ^ + 25R(e) 



\B{ 



|6P 



where Q = {AB)/{Aab). Furthermore we will rewrite |^'^(a;)p by using Jacobi polynomials. Let P^'^^{x) 
denote the Jacobi polynomial which is orthogonal on [—1, 1] with respect to (1 — xY{l +x)'^ with v^p,> — 1. 
Then the following relation holds: 



T{n + v+\) 



iFi{-n, n + v + n+l:v+l;{l- x)/2), 
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where 2^1 (a, b; c; z) is the hypergeometric series and T{z) is the gamma fmiction. In general, as for orthogonal 
polynomials, see [7]. Then we have 



\b\' 



E 

7=1 

^ \ l«|2 



2 \ 7-1 



\ (l-\\(n-l-\ 



\-2{l-\) 



I 



7=1 



2 X 7-1 



7— l/ \ 7— 1 



|„|-2(/-l)p0^«-2/(2|a|2 - 1). 



(4.4) 
(4.5) 



By using Eqs. and (|i?5|) . we see that for 1 < ? < [n/2] 



23?(e)|(p°) 



where = -P/l" ^'(2|ap — 1) (i = 0, 1). By a similar argument in [T71I18) . we obtain 



lim W^lW^^ hm I*n(^)r = !^ / 



2 /-I"! /i^(x) 



2^ 



(4.6) 



where 



|aP|6|2 

and 3ff(z) is the real part of z G C. Indeed, if n 00 with kjn e (—(1 — |a|)/2, (1 + |a|)/2), then 



pO 



2|a|2'= 


— n 






2|a|2'= 


— n 







cos(yln + P), 



■ cos(An + B 



where A = (1— |ap){(2x— 1)^ — |ap}, A and B are some constants which are independent of n, and 9 G [0, 7r/2] 
is determined by cos 6 — {\ ~ |ap)/4x(l — x). By these asymptotics and the Riemann-Lebesgue lemma, 
we have Eq. (14. 6p . 
Noting that 



\a\ 



dx 



"I l-x dx 



\a\ 



\a\ 



dx 



1 + x 



x^ 

(2-yT-H^)7r _ {2-\b\)7T 



we have 



(2-|6|)|A|2 ^ |6||P|2 



lim 



+ 



\b\ 



+ |6p(|6|-l)5R(e). 



2 2 

Moreover, the following relations hold: 

\A\' = \a\M' + \bm^+r, \B\' = ma\' + \am^ 

iH'-\b?)r 



^ie) = \/3\'-\a[ 
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where T = aabl3 + aab(3. Therefore we obtain the deshed conclusion. Concerning ||^^||, we get the following 
result by a similar fashion: 

hm pf^^^ H = i-i- / dx, 



where 



|a|2|6|2 l-x t |a|2 
So the desired conclusion is obtained. Finally we consider (5'^, ^'^) case. For 1 < ^ < [n/2] 



7=1 5=1 



1+S 



; - 1\ - 1\ /n - ^ - 1\ /n - / - 1 



l{n-l)AB I ( AB \A\^\ n - 1 f AB \B\^\ AB AB \B\^ - \A\ 



7(5Aa2 7 \Aa2 a6 



Aa2 ab J Aa2 aF 



a6 



Hence we have 



where 



lim ^-^^ ^ ' ^ ^ 



2vr 7-|a| 



dx, 



|A|2\ 1 



Aa2|6|2 VAa2 a6 / 1 + i 



AB |B|2\ 1 AB AB |B|2-|A|2 



Aa2 a6 / 1-x Aa2 A6 



a5 



This proves the case of (4',^, 4'.^). So the proof of Theorem 13. II is complete. 

5 Results on the Shannon entropy 

In this section, we compute the Shannon entropy of the QW. Let Pn{x) = P{Xn = x). Then we define the 
Shannon entropy of the QW by 

n 

Sn^-Y. Pn{x)\0g^Pn{x). 
x— — n 

Furthermore, we let 



E 



|^-D||2 



I* 



D\\2 



for D G {L,R}. In this section, we present limit theorems of the Shannon entropy for S^,S^, and Sn- 
Define 



n 



lim 



D\\2 



(5.7) 



and 



\bl_ h^jx) 
27r ^|a|2 - a;2' 
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with D e {L,R}. Moreover 



1^1 /, A„,,2 ,«,2 , r 

7r(l-a;2)^|a|2-a;2 



/(x) = r (.) + r{x) = — '^^ , 1 - |a, , 



The fohowing long-time behavior of the Shannon entropy is shown. 
Theorem 5.1 // the QW is determined by U with abed ^ 0, then we have 



lim , y = hm - — = 1, (5.8) 

n^oo log2(n/2) n^oo log2(7T./2) 



where D e {L,R}. Furthermore, 



hm^ log2 (n/2) ( ^33^^ - 1 ) = - / log2 ) ^2:, (5.9) 



hm log2(n/2) ( ^" -l]= - r f{^) log2 (/(a;)) dx, (5.10) 



w/iere D £ {L, R}. 



Bracken et ah [TO] numericahy computed the Shannon entropy of the QW defined by U with a = —d = 
yfp, b = c = 1/1 — p {p = 0.5 case is the Hadamard walk). The numerical result of Chandrashekar et al. \T2\ 
suggests that the Shannon entropy Sn of the QW given by U{9) {0 < 9 < n/2) is maximum for 9 = 7r/4 (the 
Hadamard walk), and Sn{0i) < Sn{92) < 9i < 02 < 7r/4 and Sn{9i) > Sn{92) if 7r/4 < 9i < 02 < 7r/2 
for any n. Our result is consistent with the numerical results. Moreover, we give the limiting values of the 
Shannon entropies explicitly which are new findings of this article. 

When Pn{x) is a uniform distribution on Wn — {— —{n — 2), . . . ,n}, i.e., Pn{x) — l/{n+ 1) for any 
X € Wn, then Sn takes the maximum value 

5-- = l0g2(n + l). 

In Eq. we can replace log2(ri,/2) by S'™'"'", since 5'^''''' = log2(n + 1). 

For the symmetric random walk whose walker moves one unit to the right with probability 1/2 and to 
the left with probability 1/2, the following result of the Shannon entropy, Sn^ , is known, (see [HIUS] for 
more detailed information): 

oRW 

hm , " ^. = 1. 

n->-oo log2(v^j 

Furthermore, 

/ oRW \ 1 

lim log2(v^) ■ " -1=9 (27re) . 

n^oo Vl0g2(V«) / 2 

This shows that the rate of convergence of is related to that of walker's position in the central limit 

theorem (CLT). Eqs. (j5.8p . (|5.9p and ()5.10p shows that in the QW case, the rate of convergence of the 
Shannon entropy is also related to that of the CLT. But the relativity is slightly different from the random 
walk case. 

Proof. We here consider only case. Other cases, i.e., S'^ and Sn, can also be shown in a similar way. Sj^ 
is rewritten as 

1 

^„''=log2(ll*^lP)-^iqp E log2(l*^(^)P)- 
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By Eq. ((5J)) . we have 



1 

lim 5^ = log2 (p^) - — lim ^ \^l^{x)\' log2 . 



a:— — n 



As in a similar proof of Theorem 13. 1[ we obtain 



n 



^ |VI/^(X)P l0g2(|vl/^(x)n 



x——n 




-\a\ 



/^(x)log2 (2/^(a;)/n) da: 




where a(n) ~ b{n) means a{n)/b{n) — 1 as n — >■ oo. The third equahty comes from Eq. (|4.6p . Thus 



and this completes the proof. 

6 Summary 

In this paper, we show hmit theorems for the von Neumann entropy and the Shannon entropy of the QWs 
on Z starting from the origin with arbitrary coin and initial state. In order to compute the entropies, we use 
a path counting method which is a powerful tool for analysis of QWs on Z. By using this method, we can 
derive limiting values of the entropies explicitly. But in general, it is hard to use the path counting method 
for higher dimensional cases. For the RW case, the Shannon entropy plays an important role in the theory 
of large deviation principle. It is a natural question what is a role of the entropy in the large deviation 
principle for the QW. 
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